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Abstract. We investigate the ground-state energy of the n~p atom (pionic hydrogen) in the framework of 
QCD+QED. In particular, we evaluate the strong energy- level shift. We perform the calculation at next- 
to-leading order in the low-energy expansion in the framework of the relevant effective field theory. The 
result provides a relation between the strong energy shift and the pion-nucleon S'-wave scattering lengths 
- evaluated in pure QCD - at next -to-leading order in isospin breaking and in the low-energy expansion. 
We compare our result with available model calculations. 
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1 Introduction 

The theory of strong interactions has entered a high pre- 
cision phase in recent years, both experimentally and the- 
oretically. On the experimental side, we mention i) the 
muon anomalous magnetic moment measurement at 
Brookhavcn A calculation of (g — 2) M that matches the 
foreseen experimental accuracy requires that the cross sec- 
> tt + it~ in the low-energy region is known to 



better than one percent; ii) experiments that aim to deter- 
mine hadronic scattering lengths with high precision are 
presently running at CERN [g and at PSI j|,§|. On the 
theory side, the tttt scattering lengths have recently been 
calculated at the few-percent level in Ref . || . In this arti- 
cle, we are concerned with the ongoing experiment on the 
measurement of the energy levels and decays of the ir~p 
atom (pionic hydrogen) at PSI |j,pj . It is planned to mea- 
sure the strong interaction width and shift of the ground 
state at the percent level. These measurements can then 
be used to directly extract from data the ttN scattering 
amplitude at threshold. The aim of the experiment goes, 
however, further: it intends to extract from these measure- 
ments the S'-wave ir N scattering lengths + al + in pure 
QCD with high precision. In order to achieve this goal, the 
relation between the scattering lengths and the threshold 
amplitude must be known at an accuracy that matches 
the accuracy of the experiment. In other words, one has to 
remove isospin-breaking effects from the threshold ampli- 
tude with high precision, in the framework of QCD+QED. 



Isospin-breaking effects in low-energy wN scattering 
have been extensively discussed in the literature on the ba- 
sis of a potential model approach. The discussion relevant 
to the problem of pionic hydrogen can be found in Refs. j^, 
0| . Here, we rely on the effective theory of QCD+QED, a 
method already invoked p,p[ |ic| , pl| in the analogous pro- 
gram for the Tr + n~ atom investigated presently at CERN 
For a critical comparison between the potential model 
and the effective theory framework, we refer the reader 
to Ref. fl. 
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In the case of pionic hydrogen considered here, the 
relation between the strong energy-level shift e\ s of the 
ground state and the threshold n~p scattering amplitude 
has been worked out in the effective theory in Ref. , 

€ls ~ 2ttM„ \ 2nM„ J ' ' ' ' 

(1.1) 

Here, /i c = m p M 7T (m p + M^)' 1 denotes the reduced massQ 
of the tt~p system, and T^n is the threshold amplitude 
for the process ir~p — > n~p, evaluated at next-to-leading 
order in isospin breaking. Furthermore, a ~ 1/137.036 
denotes the fine-structure constant, and the ellipsis stands 
for the higher-order terms in isospin breaking (see below). 
In the isospin symmetry limit, the threshold amplitude 
is proportional to a particular combination of the S'-wave 
pion-nucleon scattering lengths (we use the notation 
of Ref. @), 

% N = T" N + aT~< + (m d - m u )T m , 

T° N =4^1+ ¥A (a+ + + a - + ) . (1.2) 
V m p / 

The isospin breaking amplitudes T 7 ' m depend on the 
renormalization group invariant scale of QCD, and on 
the quark mass rh — (m u + md)j1 (we consider two- 
flavor QCD). Once T 7,m are calculated, a measurement 
of the shift ei s allows one to determine the combination 
Oq + + aw . Similarly, a width measurement of the ground 
state provides \% + \. The values of these scattering lengths 
are correlated e.g. with the pion-nucleon sigma-term Jig ], 
with the pion-nucleon coupling constant fll6| , and with 
the Goldberger-Treiman relation, which relates the pion- 
nucleon coupling constant to the axial charge. It goes with- 
out saying that the scattering lengths are therefore very 
central objects in the analysis of pion-nucleon reactions. 
We refer the reader to M for a recent investigation of these 
questions. Here, we concentrate on the evaluation of the 
isospin breaking amplitudes T 7,m . 

The amplitudes T 7 ' m can e.g. be evaluated in the 
framework of chiral perturbation theory. Writing the en- 
ergy shift in the form 

e ia = -2a 3 f i 2 c (a+ + + o^)(l + S e ) , (1.3) 

the leading order calculation |l| gives S e = (—4.8 ± 2.0) • 
10 -2 , whereas the calculation by Ref. Q in a potential 
model framework leads to 8 e = (—2.1 ± 0.5) • 10~ 2 . How 
should the two calculations be compared? As has been 
pointed out in Ref. the leading order terms in the 
effective theory are due to effects that are not all con- 
sistently taken into account in the potential model cal- 
culation. On the other hand, as has been emphasized in 
footnote 1 in Ref. j| , mass splitting effects in strong loops 
and jn intermediate states show up only at higher orders 

1 We denote the charged pion mass (the proton mass) with 
Mi(m p ). 



in the chiral expansion. It is the aim of the present article 
to carry out the calculation of T 7 and T m at next-to- 
leading order, where these effects come into play. Further, 
the interference effect between vacuum polarization and 
strong interactions will be taken into account as well|^] - 
they occur formally at next-to-next-to leading order in 
isospin breaking. 

Isospin breaking effects in ttN scattering have been 
studied already previously [see Ref. Q and references 
therein] by using heavy baryon chiral perturbation theory 
(HBChPT). We were not able to directly use these results, 
for the following reason. In Ref. p^ , the authors calcu- 
late physical amplitudes in different channels, and study 
combinations thereof that vanish in the isospin symme- 
try limit. This is not what is needed for pionic hydrogen, 
where one has to extract isospin-breaking correction to a 
single amplitude (tt~ p e lastic amplitude, in the case of the 
energy levels). Ref. |l7| does not provide explicit expres- 
sions for the physical amplitudes, and we have therefore 
performed an independent calculation. 

In Ref. |l7| , HBChPT was used to calculate the am- 
plitudes. Here, we rely on the framework developed by 
Becher and Leutwyler [18|. This method is manifestly 
Lorentz invariant and preserves chiral power counting in 
the case where baryons and Goldstone bosons are running 
in the loops. In the present context, contributions from vir- 
tual photons are needed as well. We show that the method 
proposed in Ref. Q can straightforwardly be adapted to 
this case. Individual Feynman diagrams contain ultravi- 
olet (infrared) singularities, due to integration over large 
(small) momenta. We use dimensional regularization to 
tame both of these singularities. The counterterms from 
the higher order chiral Lagrangians cancel the ultravio- 
let poles at d = 4 in the final result. In order to check 
this cancellation, we evaluate the 1-loop divergences also 
with the heat-kernel method. Finally, infrared divergences 
disappear in physical quantities at the end of the calcula- 
tion. In our case, this concerns the n~p — ► 7r~p scattering 
amplitude at threshold, with the Coulomb singularity re- 
moved. Lack of phase space forbids the emission of soft 
photons, as a result of which the elastic scattering ampli- 
tude must be infrared finite at threshold. Needless to say 
that this cancellation serves as another welcome check on 
our calculation. 

The final result for S e contains, at the next-to-leading 
order considered here, several low-energy constants 
(LECs) that parameterize the structure of the effective 
theory at this order. All but one of these LECs are un- 
der experimental control or expected to generate a small 
effect. The remaining one, /i, enters the expression al- 
ready at leading order [|l3). Whereas this constant is ex- 
pected to contribute sizeably to S e , no precise determi- 
nation in terms of experimental data is available yet. We 
shall explain why, in our opinion, a precise determination 
of aj + + (Iq + from a measurement of the energy-level shift 
t\ s has to await a corresponding precise determination of 



2 We are indebted to T. Ericson for pointing out to us that 
this may be an important effect. 
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/i, and more importantly, why potential models are not 
of help in this respect. 

The layout of the paper is as follows. In section || we 
introduce the definition of the threshold scattering ampli- 
tude in the presence of photons, and give the relation of 
this quantity to the strong energy shift of the ir~p atom. 
In section we display the set of the meson and meson- 
nucleon chiral Lagrangians, which are used in the calcu- 
lations of the threshold scattering amplitude at 0(p 3 ). In 
section | we calculate the tree-level contributions to the 
threshold amplitude at <3(p 3 ). In sections |]and^, we out- 
line the generalization of the infrared regularization pro- 
cedure Q, needed in the presence of virtual photons. The 
general procedure for the calculation of the ^-matrix el- 
ements in the infrared regularization in the presence of 
virtual photons is described in section |?], where we also 
discuss the separation of infrared and ultraviolet diver- 
gences. Section || deals with the calculation of the vertex 
diagram with virtual photons. Here, we demonstrate the 
origin of the Coulomb phase and the singular piece of the 
scattering amplitude at threshold. In section ^ we present 
the final result of our calculations of the threshold am- 
plitude, and in section ^ we discuss the size of the rele- 
vant low-energy constants. In section [H], we evaluate the 
ground-state energy-level shift numerically and compare 
the result with model calculations. Section |l2] contains our 
summary and conclusions. Technical details are relegated 
to the appendices. 



2 Threshold amplitude and the strong 
energy-level shift 

It is both conventional and convenient to split the ground- 
state energy into electromagnetic and strong parts, 



Eis — + ei 8 



(2.1) 



Despite the fact that this splitting cannot be understood 
literally (e.g. there are contributions from strong interac- 
tions in Ei s n ), it turns out to be useful in the theoretical 
analysis of hadronic atom data. 

The electromagnetic part has been worked out in 
Refs. The results of these two investigations - that 

were performed in completely different settings - agree nu- 
merically with high precision. The complete analytic result 
of the electromagnetic energy at 0(a 4 ), which is given in 
Ref. jlj|, includes relativistic corrections, finite-size effects 
from proton and pion charge radii, and the correction due 
to the anomalous magnetic moment of the proton. In ad- 
dition, it includes the contribution from the electron vac- 
uum polarization, which arises formally at 0(a 5 ), but is 
amplified by a large coefficient (M^ / m e ) 2 . The evaluation 
of Ef^ in Ref. || includes in addition some higher-order 
corrections (vertex corrections, corrections to the vacuum 
polarization contribution), that were omitted in Ref. pH| . 
Numerically, however, they are much smaller than the ones 
included in both approaches. 

What is exactly measured in the experiment at PSI, 
is the transition energy E^^l s between the 3p and Is 



states Since the electromagnetic shifts of both lev- 

els are known to a very high precision and since the 
strong shift in the 3p state is very small, the measurement 
allows one to directly extract the strong shift of the Is 
level according to 



iriem 



pmcas 
-"3?- Is 



(2.2) 



The aim of the new experiment at PSI is to measure e\ a at 
the percent level Once this measurement is performed, 
one is faced with the task of extracting the strong S*-wave 



wN scattering lengths a 



- t o+ 



from this quantity. 



In order to formulate the problem rigorously, it is 
convenient to introduce a common counting rule for the 
isospin-breaking effects, parameterized by the fine struc- 
ture constant a and the up and down quark mass differ- 
ence md — m u . We have found it useful to count these ef- 
fects at the same order, intr oduc ing the formal parameter 
S ~ a ~ md — m u . Equation (1.1), which relates the strong 
shift e\s to the threshold amplitude, is then valid jl3| up to 
and including terms of order S 4 - hereafter, this is referred 
to as the "next-to-leading order in isospin breaking" . In 
the following, it is useful to introduce the notation 



%tN — T® N + ST , 
ST = aT 1 + (m d - m u )T r ' 



(2.3) 



The isospin-symmetric part T® N refers to pure QCD, 
where the pion and the nucleon masses are equal, by con- 
vention, to the charge d pio n and to th e pr oton masses, re- 
spectively. From Eqs. (|l.l[), ( |l.2| ) and (1.3), the expression 
for the isospin-breaking correction S e can now be readily 
worked out. One still has to add the correction <5J ac due 
to the interference of vacuum polarization and strong in- 
teractions. The complete expression takes the form 



ST 



47r(l + M 7r /m p )(a^ + + a 0+ ) 
K = — 2a(lna — l)/z c (aj , + ag+) 



K + 5™ c , 



(2.4) 



The only unknown ingredient in Eq. (2.4) is thus the 
isospin-breaking part ST of the threshold amplitude. The 
chiral expansion of the amplitudes T 7 and T m is per- 
formed in the variable r = M^/m p - the expansion of ST 
therefore starts at order p 2 (counting a and md — m u as 
0(p 2 ), as usual), 



ST = ST 2 + ST 3 + 0(p 4 ) . 



(2.5) 



The leading term ST has been determined in Ref. [ p3| . 
Here, we evaluate 5T 3 , which amounts to a one-loop cal- 
culation of the n~p — > n~p amplitude. We will show be- 
low, that the chiral expansion of T m starts at 0(r 2 ). At 
the accuracy we are concerned with here, the quark mass 
difference does thus not enter. [The effect coming from the 
insertion of the quark mass difference into loops may be 
non-negligible E9|. Its evaluation requires the calculation 
of the threshold amplitude at order p 4 , which is beyond 
the scope of this work.] 
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The contribution <5J ac has been omitted in Ref. fll3[ , 
since formally it is higher order in the parameter 6. How- 
ever, numerically it is not negligible, because of the large 
coefficient containing M^/rrie. The calculation of this cor- 
rection within a non-relativistic effective Lagrangian ap- 
proach has been carried out in Ref. [20[. Here, we use this 
result, contained in table II of Ref. 





tf(0) 



= 0.48 % 



(2.6) 



where ^(0) stands for the pionic hydrogen wave function 
at the origin, and where 6\P(0) denotes the correction to 
this wav e function due to vacuum polarization. The re- 
sult (|2.6| ) agrees perfectly with the value (5J ac = 0.46 % 
obtained within the potential model of Ref. || . 

We now present the exact definition of t he t hreshold 
amplitude T^n which enters the expression (1.1) for the 
strong energy shift of the ground state. Let us consider 
the elastic scattering process 7r~(g) + p{p) — » n~(q r ) + 
p{p') in the vicinity of the physical threshold, at first order 
in the fine structure constant a. External momenta are 
on the mass shell p 2 = p' 2 = nip, q 2 = q' 2 = M 2 . The 
Mandelstam variables are defined in the standard manner, 
S = (p + q) 2 = (p' + q') 2 , 
U= (p - q') 2 = (p> - q) 2 , 

The 3-momentum of the pion and of the nucleon in the 
CM frame is given by |p| = A 1 / 2 (s, m 2 , M 2 )/(2^/s), where 
X(x, y, z) — {x—y— z) 2 — Ayz denotes the triangle function. 
The scattering angle is given by cos6> = 1 + t/(2p 2 ). 

The scattering amplitude for this process can be ex- 
pressed in terms of two scalar functions D and B^ 



t = (P - v'Y = (? ~ q') 2 , 
s + t + u = 2(m 2 p + Ml). 



T 



u(p'){D(s,t) 



4m, 



[4\$]B{s,t)}u{p) 



u(p)u(p) = 2m f 



(2.7) 



In the context of the energy-shift considered here, only the 
amplitude D(s,t) is relevant. In the next step, we define 
the truncated amplit ude Z), which is obtained from the 
scattering amplitude (2/7) by subtracting the one-photon 
exchange contribution displayed in Fig. [y, 



D(|p|,cos0) = L>(s,i) = D(s,t) 



2 ^(f)F!(t)( S - U ) 



2m p t 



(2.8) 



where F w (t) and Fi(i) denote the pion electromagnetic 
form factor and the nucleon Dirac form factor, respec- 
tively. 



3 It is well known that the above scattering amplitude in the 
presence of electromagnetic interactions is infrared-singular in 
perturbation theory. In this paper we use dimensional regular- 
ization to tame both, infrared (IR) and ultraviolet (UV) diver- 
gences. Thus, the scattering amplitudes D and B in Eq. (2.7) 
are meant to be evaluated at d ^ 4 (see below). We use the 
notation of Bjorken-Drell |2l] for the Dirac matrices. 



Fig. 1. Definition of the truncated 7r~p — + 7r~p amplitude. 
The filled vertexes in the diagram with the photon exchange 
correspond to vector form factors of the pion and of the proton, 
calculated in pure QCD. T w m corresponds to the truncated 
amplitude. 



In the presence of virtual photons, at order a, the scat- 
tering amplitude -D(|p|, cos 9) contains Coulomb singular- 
ities at the physical threshold |l3|,[ll]]. In particular, there 
emerges the Coulomb phase which is divergent in physical 
dimensions d = 4 (cf with p2[), 

- I(r'(l)+ln47r)+ln^Ml . (2.9) 
2 fx J 

Here fi stands for the scale in dimensional regularization. 

After removing the Coulomb phase, the behavior of 
the real part of the amplitude Z?(|p| , cos 9) in the vicinity 
of threshold, at order a, is given by |13| 



Re{( 

P 



-2ia0c ( |p| 



D(|p|,cos0)} 



|p-o 



B, In ■ 



He 



N 



0(\p\), 



(2.10) 



where the quantities B\ , Bi can be related to the S'-wave 
ttN scattering lengths (their explicit form is not needed 
in the applications to the hadronic atom problem). The 
quantity 7^, which is referred to in this article as the 
threshold scattering amplitu de, is infrared finite. It deter- 
mines ei s according to Eq. (1.1), and is the central object 
of our investigation. 



3 Lagrangians 

In this section, we display the low-energy effective La- 
grangians which are used in the calculation of the ir~p 
amplitude. In particular, we provide the meson-nucleon 
Lagrangian at order e 2 p, whose divergences are consis- 
tent with the "standard" choice of the meson 0(p 4 ) | p3[ , 
0(e 2 p 2 ) [Q, and nucleon 0(p 3 ) p5| , p6| Lagrangians. The 
reason for doing so is that the 07/e~pT Lagrangian which 
is available in the literature [Jl^, g^,|2g] , is not consistent 
with this choice. The details of the derivation, which was 
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performed by using the Berezinian approach [ [29| , |3C[ ] , can 
be found in appendix [b]. In this appendix, we in addi- 
tion give the relation between two sets of the LECs which 
are defined by the choice of the Q{p A ) meson Lagrangian 
either in form of Ref . |23| , or Ref . |3l[] . 

The full Lagrangian consists of meson and nucleon 
parts, as well as the free photon Lagrangian together with 
the gauge-fixing term, 

^■eff = £tt + £jV + £7 ; 

r - rip 2 ) + r(e 2 ) + rip 4 ) + r( e V) , . . . 



£ 



N 



— £^ - 1 - r< ~ p ' 



c 



TV 



(e 2 ) 
JV 



(e 2 p) 
JV 



(3.1) 



where 



F 



C<* ) + £i e } + £ 7 = — (d"fftd„ff + x f ^ + rf X ) 



-ZF\QUQU ] ) - J F tlv F ,a ' - ^{d^f , (3.2) 



(3.4) 



where U is a unitary 2x2 matrix, and \P is the nucleon 
field. As usual, s,p,v^,a^ denote external scalar, pseu- 
doscalar, vector and axial fields, B is a constant related 
to the quark condensate, and Q — ediag(|,— i) is the 
quark charge matrix. The axial source is taken traceless, 
K) = 0. 

Building blocks for the meson-nucleon Lagrangians in 
the presence of virtual photons are 

Dfj, = d,j, + r M , U = u 2 , — iu^d^U v) , 
F^ = - [u\ d^u] - % - v)R^u - l - uLpU 1 , 

x± = ti t xu t ± ux^u , x+ = x+-^(x+), 

Q± = \ (uQu f ± u^Qu) , Q± = Q± - ~ (Q±) , 



7 10 

4«V) =F 2^ fc . (eV) ) 
i=l i=l 



i=l 
23 



12 



(3.3) 



!=1 



1=1 



where (A) denotes the trace of the matrix A. The building 
blocks for the mesonic Lagrangians are 



£„ 



v t _ l ±a t ^ + QA,j. , 



= vtU^u ± uC^u 1 , 
• / "± — At _ 2 ' ± ' ' 

X = 2B(s + ip) , d^x = d^x - iR^X + *x£^ > 
Q± = -{uQu^±u^Qu), 



(3.5) 



where Q — ediag(l,0) denotes the nucleon charge ma- 
trix, and Lp, R^, L^ v are defined just like their pi- 
onic counterparts TZ^, £^, TZ^, L^ v respectively, with Q 
replaced by Q. We have set the charge matrices Q and 
Q to their constant physical values, and do not consider 
the most general expression of the effective Lagrangians 
containin g sp ace-dependent spurion fields Ql,r- Further- 
more, in (|3.3| ) we drop terms which do not contain pion 
fields, and terms of order e 4 . 

We comment on the LECs F, gA, Z,a, fi, ■ ■ -. The first 
one, F, denotes the pion decay constant in the chiral limit, 
and gA is the axial charge, again at m u — md — 0. The 
quantity Z is related to the pion mass difference, 

M 2 = (m u + m d )B + 0(m 2 q , e 2 m q , e 4 ) , 

Ar = M 2 - M 2 = 2e 2 F 2 Z + 0(m 2 q , e 2 m q , e 4 ). (3.6) 

The couplings Cj, /, are finite, if the calculation of the loop 
diagrams is performed in a manner that respects chiral 
power counting. The divergences in the remaining LECs 
are given by 



di = A + <(//), gi = j^X + glfa), (3-7) 



where 



A = 



d-4 



1 1 



16tt 2 2 



--[r'(l) + ln47r + l] . (3 



In appendix ^ we list the operators o[ k ^ , as well as the 
divergent parts of the low-energy constants, and provide 
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references to the original literature. Since our operator ba- 
sis in the mesonic sector is defined in a standard manner, 
numerical values of 1%, k\ can be directly taken from the 
existing analyses. The definition of the finite constants Cj 
and /, coincides with that from Refs. p5| and p7||, re- 
spectively (in the latter, the notation /^instead of /$ is 
used for the LECs in the relativistic theory). Our choice 
for those d\ that contribute to the irN scattering ampli- 
tude, corre sponds to the one of Ref. modulo the rela- 
tions (Bf4), which display the effect of the different choice 
of the mesonic basis here and in Ref. ]3^] . Finally, for the 
reasons given in appendix |b| it is not clear to us how we 
can compare the couplings g\ used here, with the ones 
determined in Ref. 



;oupl 

m 



In the actual calculation of the scattering amplitude, 
we invoke an exponential parameterization of the matrix 
U, 



U = cupiir/F , 7r 



-( 



\V2tt- 



V / 2tt + 



(3.9) 



and use 7r as an interpolating field for the pions. Further, 
we set 



X = 2Bdiag(m u ,m rf ) , 



(3.10) 





Fig. 3. Counterterm contributions. Insertions on the external 
lines are not shown. 



The scattering matrix elements are evaluated with me- 



12 



Mi 



son momenta that are on the mass shell, q 2 = q 
and with spinors that obey the Dirac equation with the 
physical proton mass, jju(jp, s) — m p u(p, s). On the other 
hand, the neutron mass in the internal line in Fig ||b is the 
chiral symmetric one. It receives corrections from coun- 
terterm insertions © in Fig^ 3|, and from the self energy 
diagrams si and S2 in Fig. O. Here, we merge these con- 
tributions and evaluate the diagram Fig ^b with the phys- 
ical neutron mass. [The subtraction in the one-loop self- 
energy-type diagrams si and S2 in Fig. [ll] is then carried 
out on the mass shell.] In this convention, the contribution 
from the Born diagrams Fig.0 becomes 



T = 



2F 2 



2F 2 (r 



(4.1) 



and perform the calculations in the Feynman gauge a = 1 . 



4 Tree contributions 



The tree contributions to the threshold scattering ampli- 
tude n include the pseudovector Born diagrams Fig. p| 
and the contributions from the counterterms in the ef- 
fective Lagrangians, indicated in Fig |3[ Insertions in the 
external lines are not displayed - they will be included 
through mass and wave function renormalization in sec- 
tion pL 



Next consider the vertex corrections © - these can be ab- 
sorbed in g^. Since they are proportional to the quark 
masses m u ^m^ they do not contribute at order p 3 , and 
we drop them altogether. The polynomial contributions© 
arise from two-nucleon two-pion vertexes in the effective 
Lagrangians. We finoQ 

T« = -i^e 1 + M (c2 + Ca) _^ (4/l+/3) 



4M 3 

(di+d 2 + 
+2e 2 M 7r (g 6 +g 8 ) 



d 3 ) 



8AUM 2 



F 2 



^d 5 



5 Two-point functions of pion fields 



(4.2) 



(*>) 



Fig. 2. The Born diagrams, vector (a) and axialvector (b) 
couplings. 



In Ref. |Tq ] , Becher and Leutwyler have shown how power 
counting in baryon chiral perturbation theory (pions and 
nucleons running in the loops) can be incorporated in 
a manifestly Lorentz invariant manner. Here, we extend 
this framework to include virtual photons. In addition to 
the question of power counting, virtual photons generate 
the standard problems: poles in two-point functions are 
transformed into branch points, wave function renormal- 
ization constants become ill-defined, and truncated on- 
shell Green functions in general cease to exist in four 



4 The normalization of the low-energy constants /^j used in 
the present paper differs from that of Ref. O: F 2 f^ c = f?$. 



J. Gasser et al.: Ground-state energy of pionic hydrogen to one loop 



7 



dimensions. In order to identify these infrared singular- 
ities 5 , we start the discussion of the low-energy expansion 
of Green functions with the two-point function of pseu- 
doscalar quark currents. In this case, power counting does 
not pose a problem, because one is concerned with mesons 
and photons only. Furthermore, ultraviolet divergences do 
not show up in the final result either, because Green func- 
tions of quark currents are well defined objects in the ef- 
fective Lagrangian framework. Therefore, the only obsta- 
cles in this case are the infrared divergences generated by 
photon loops. 



5.1 Pseudoscalar two-point function 




Fig. 4. Photon contribution to the two-point function of pseu- 
doscalar fields. 



The pseudoscalar densities P^(x) — q(x)ij5T 1±l2 q(x) 
may be used as interpolating fields for the charged pions. 
To evaluate scattering matrix elements, the residue of the 
two-point function 



G(s)=i / ^e-^(0|TP 5 -( a; )P 5 + (0)|0) 



s = q 



(5.1) 

at s = M 2 is needed - it determines the wave function 
renormalization constant. However, in the presence of vir- 
tual photons, G develops a branchpoint, not a pole, and 
that residue does not exist. There are two standard proce- 
dures to deal with the situation: Either, one introduces a 
photon mass, which shifts the branchpoint to s — (M n + 
m 7 ) 2 , as a result of which G indeed contains a pole at 
the pion mass. The photon mass is then sent to zero at 
the very end of the relevant calculations of physical quan- 
tities. Or one instead uses dimensional regularization to 
tame both, ultraviolet and infrared divergences. In the 
following, we use the latter method. Let us consider the 
low energy expansion of G. Including contributions from 
one virtual photon as indicated in Fig. 0, and adding the 



/ 2 2 \ 

counterterms from p in Eq. (3.1), we find 



G(s) 



AF 2 B 2 R{s) 
AB-s 



16tt 2 



2 Jl/r2 j^r 
M 



e M K 



M 2 = (m u + m d )B , 



M 2 
V 2 



is) 



s + M 2 . M 2 



In- 



M 2 



2f) 92 4 

K r — <h r _l_ h r h r \ — —h r — —h r 



v 8 i 



There are two types of infrared singularities in the present 
context: The ones associated with non-analytic terms in the 
chiral expansion, and singularities generated by the presence 
of photons. Since it will always be clear what singularities we 
have in mind, we do not distinguish in the following between 
the two, which makes notation less clumsy. 



on o 

K r G = T (fcl + k r 2 - 2kl - 2kl) - -k r 7 - 8k r a . (5.2) 
[Here, we have omitted the contributions from the lead- 



ing electromagnetic term proportional to Z in Eq. (3.6), 
such that M 2 = M 2 at leading order.] The ellipses de- 
note additional terms in the low-energy expansion. The 
correlator G(s) is finite at nonexceptional momenta, but 
the residue develops a logarithmic singularity at s = M 2 , 
manifest in the function r(s) [at the order of the expan- 
sion considered here, we may replace M 2 in the function 
r(s) by M 2 }. 

In order to continue, we keep the dimension different 
from four until the very end of the calculation of physical 
quantities. The function r(s) becomes 

r( a ,d) = < ?-^(Aiiyr{uj)M- 2 " [ dxx-" 



o 



a - 1 

x [ X - w -(1 - a + ax)'"] + 0(d - 4) 
oj = 2 - - , a = s/M 2 . 



(5.3) 



At nonexceptional momenta, r(s, d) approaches r(s) as 
d — > 4. If d is bigger than four, one may perform the 
mass-shell limit, 



R(s,d) = l-^ 



— C 



ID 



In 



M 2 



e 2 K r G + 0{d - 4) 



Mi 



d-4 



--(r'(l)+ln(47r) + l) 



(5.4) 



The infrared singularity manifests itself in a pole at d = 4, 
which is due to integration over small loop momenta. We 
indicate the origin of this singularity with the symbol IR 
in the divergent quantity C^. 

We do not use pseudoscalar densities as interpolating 
fields - it is simpler to use the pion fields instead. In this 
case, additional singularities occur at d = 4, as is shown 
in the following subsection. 
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5.2 Two— point function of the pion fields 



The Green functions of the pion fields depend on the pa- 
rameterization used for the matrix U(x) in the effective 
Lagra ngia n. Here, we use the exponential parameteriza- 
tion ( |3.9D , and consider the propagator of the charged 
fields, 



R-*(s,d) 
Ml - s 



dxe 



|Tvr-(x)7r + (0)|0) 



s = q 



(5.5) 



We again include the photon loop from Fig. and the 

pertinent counterterms from P at Z = 0. The result 
for the pion mass is the same as before, whereas several 
of the counterterms fc, are now absent in the residue. We 
find 



R-n(s, d) = 1 — 



32tt 2 



M 2 

4r(s,d) + C^ W + 31n— -4 



uv 



2M. 



d-4 



10 

Kl = -{k\ + kl) , 
(^4 ^(m + ln(4.) + l)) 



(5.6) 



In contrast to the two-point function G(s), the propagator 
of the pion field is ultraviolet divergent: the function R v 
develops a pole in four dimensions also at nonexceptional 
momenta. We indicate the origin of this singularity with 
the symbol UV in the divergent quantity Ctjv- The rea- 
son for the occurrence of this pole is well understood: the 
effective theory guarantees that Green functions of quark 
currents are ultraviolet finite, whereas the pion fields sim- 
ply serve as integration variables in the path integral and 
are devoid of any physical significance. 

We conclude that dimensional regularization generates 
two different types of poles at d = 4: those that are due to 
integrations over large loop momenta, and those that are 
due to considering the mass shell restriction (integrations 
over small loop momenta). In each graph, one may easily 
distinguish between the two types of singularities, and we 
will do so below, by denoting the corresponding singular 
quantities with indices UV and IR, respectively. As we 
shall see later, distinguishing between the two divergences 
serves as a powerful check on our calculation. 

Later in this article, we also need the contributions 
from the pion loops displayed in Fig. 5. Including also 
the terms generated by the leading electromagnetic term 
proportional to Z which contributes to the pion mass dif- 
ference, we find 

Z, = RAMld) = 1 + -fl (C^ - d v ) 



+ 



1 



24 7T 2 F 2 



{Ml 



8tt 2 



z 2 F 2 



Z)Cl 



uv 



24tt 2 



20 e 2 



(fci + k 2 ) + 0(e 2 M 2 , M 4 , e 4 ) + 0(d - 4) 



Fie. 5. Pion loop contributions to the two-point function 
(5.5). Neutral and charged pions are running in the loop. 



(5.7) 

Here, it is convenient to use the unrenormalized couplings 

&1.2- 



6 Two-point function of baryon fields 

6.1 Self-energy of the heavy scalar field 

We now consider the case where a heavy particle interacts 
with the photon. To illustrate the method to preserve chi- 
ral power counting in this case, we consider the self energy 
diagram Fig. |(| which represents a heavy scalar particle 



x 



J. 



Fig. 6. Self-energy diagram for a heavy field of mass m. The 
solid line denotes the heavy field, and the dash-dotted line 
stands for a massless scalar particle. 



of mass m, interacting with a massless field, 
1 f d d k 



S{s,d) 



i(2w) d J [~k 2 ][m 2 - (p- kf 

= (47r)- (2 - lu) r(w) f dxz' w , 
Jo 

z = x(m 2 — s(l — xj) ; s = p 2 . 



(6.1) 



Evaluation in the standard manner shows that S is of chi- 
ral order zero, 2J ~ const., near the mass shell, whereas 
power counting suggests that it is of order one, S ~ 
(m 2 — s)/m 2 . In Ref. 18 , it has been shown that power 
counting with nucleon and pion fields is preserved, if one 
extends the integration over the Feynman parameter - that 
combines the heavy and light fields - from zero to infinity. 
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Proceeding in the same manner for photon fields, we find 
that 



S : (s,d) = (4tt)- 



(2- 



>r(u;) 



dxz 



[^)-{^)r( w )Q 1 - 2w s- 



£} = (m — s)/s . 



dxx~ w (l+x)- 



(6.2) 



where the index I denotes the infrared part a la Becher- 
Leutwyler. It is seen that the self energy diagram is now of 
chiral order one near threshold. The mass shell constraint 
is performed by first continuing the result to dimensions 
bigger than four, as before. One finds that 

£ 7 (m 2 , d) = ^'{m 2 , d) = , d > 4 , (6.3) 

where the prime denotes a derivative with respect to s. In 
other words, the self energy diagram Fig.||does contribute 
neither to the mass nor to the residue of the two-point 
function in this prescription. 

We found it very useful to distinguish between in- 
frared and ultraviolet divergences also here. This may be 
achieved as follows. We conside r th e self energy and its 
derivative at s = \ 



From Eq. (6.1), one has 



S s (m 2 ,d) = — -(-CW + 1) 

167T Z 



Ck = 2m d - i (-^ -\ (r'(l) + ln(47r) + 1)) 



K = UV, IR . 



(6.4) 



We have indicated with the index S the standard proce- 
dure, where the integral over the Feynman parameter is 
performed from zero to one. We have booked the singu- 
larity at d — 4 as an ultraviolet one for obvious reasons. 
The derivative is 

Z S \™\d) = —^ I {C m -l) , (6.5) 

where the singularity at d = 4 is now due to integra- 
tions over small momenta and therefore booked as an in- 
frared singularity. The Becher-Leutwyler infrared regular- 
ized part is obtained by subtracting from these expressions 
the 'regular part' of the Feynman diagram, 



S R (s,d) = -(4ir)- {2 - w )r(w) 



dxz 



(6.6) 



This integral converges for d < 3 - the relevant poles at 
d = 4 are therefore booked as ultraviolet. We find 



S R (m 2 ,d) -- 
S R '(m 2 ,d) 



1 



r(-Cuv + l) 



1 



32n 2 m 2 



uv 



1) 



and therefore 

S\m 2 ,d) 

S T '(m 2 ,d) 



S s (m 2 ,d) 
1 



327r 2 m 2 



S R (m 2 ,d) 
(Cm — Cuv) ■ 



0, 



(6.7) 

(6.8) 
(6.9) 



The result for the val ue of the self energy on the mass shell 
agrees with Eq. (6.3), whereas its derivative becomes the 
same upon identifying Cir with Cuv ■ 



6.2 The nucleon propagator 

The infrared regularized nucleon propagator is evaluated 
in a completely analogous manner, applying infrared reg- 
ularization to both, the pion and the photon field. In the 
following, we only need the wave function renormalization 
constant for the proton field. We define it in the following 
manner. The propagator is 

S af3 ( P ,d) = i [dxe tpx (0 1 TP a (x) P p (0)\0) , (6.10) 



where P{x) denotes the proton field. The Z-factor is ob- 
tained from 



< lim (m p - $)S(p,d) > u(p, r) = Z N u(p,r) ,d>4 , 

(6.11) 

where m p denotes the physical proton mass, u(p,r) is a 
spinor in d dimensions with (ft — m p )u(p,r) = 0, and 
where Zjq depends on the prescription used for performing 
the integration over the Feynman parameters (standard, 
regular or infrared). We obtain for the infrared regularized 
Z-factor at order p 2 



Zn — 1 + (Cm 



9A™ 2 I rrl 

Mir 2 F 2 1 



Cuv) 

V^K \ -Ouh. (0.12, 



with 



m 

Hi 



r 2 (6 



-(5- 



^9 Cuv 
3 Cuv ■ 



- 18 lnr) 
6 lnr) ; 



(6.13) 



Here, we have replaced the physical nucleon mass with its 
chiral limit mass m - this is correct at the order of the 
low-energy expansion we are working. We will do this in 
all loop amplitudes considered below. 




Fig. 7. Photo n and pion loop contribution to the two-point 
function (6.1C) of the proton. Neutral and charged pions are 
running in the loop. 



10 



J. Gasser et al.: Ground-state energy of pionic hydrogen to one loop 



7 Evaluation of S'-matrix elements 

Here, we describe the evaluation of the S'-matrix element 
for the process ir~p — > ir~p in this framework. The basic 
element to be evaluated first is the truncated four-point 
function 

(0\Tn-(x 1 )7r + (x 2 )P a (x 3 )P' 3 (x i )\0} (7.1) 

in d space-time dimensions. We take into account tree 
graphs evaluated with C e f / and one-loop graphs gener- 
ated by £i p ^ , C^k ' and £ y ] . The relevant diagrams are 
displayed in Figs. (|,|,|, 00) ■ At the end, one multiplies 
the result with the incoming and outgoing spinors for the 
proton field, and incorporates the contribution from the 
wave function renormalization constants for the pion and 
nucleon field, determined in sections |^, ^| in the standard 
manner. For general momenta of the incoming and out- 
going particles, the limit d — > 4 does not exist due to the 
infrared divergences which are generated by virtual pho- 
tons. As we discussed in section |^, these divergences are 
absent in the amplitude at threshold, provided that one 
has removed the infrared divergent Coulomb phase. Here, 
we are inter ested in the infrared finite quantity T^n de- 
fined in Eq. ( 2.10| ). The final result depends on the manner 
the integration over the Feynman parameters is carried 
out (standard, regular or infrared). Here, we work with 
the prescription given in Ref. [fl8| , which preserves chiral 
power counting. In order to keep track of UV and IR di- 
vergences, we found it convenient to evaluate the regular 
part of each diagram separately, and subtract it from the 
same diagrams calculated in the standard manner. The 
procedure was illustrated in sections ^, ^ for the two-point 
functions. It can straightforwardly be generalized to any 
of the one-loop diagrams encountered here. The schematic 
prescription is as follows: 

i) Graphs with no virtual nucleons need no infrared rcg- 
ularization - graphs with closed nucleon lines vanish in 
infrared regularization. 

ii) Graphs with pions (or photons or both) and nucle- 
ons running in the loop, are evaluated in the follow- 
ing manner. Heavy (nucleons) and light (pions and 
photons) lines in the diagrams are combined together 
separately by using Feynman parameterization, as de- 
scribed in section 6 of Ref. 



The pertinent denom- 



inators are then combined with a single Feynman pa- 
rameter [x attached to the heavy and (1 — x) to the 
light line]. The infrared regularization consists JT8| in 
extending the region of integration over this parame- 
ter, from [0, 1] to [0, oo]. 
iii) In order to distinguish between ultraviolet and infrared 
divergences in graphs with a virtual photon, we evalu- 
ate the infrared regularized part as the difference be- 
tween the standard and regular contribution. Symbol- 
ically, 



Here, M 1 stands for the final result of the calcula- 
tion of any diagram in the infrared regularization, and 
F denotes the corresponding integrand, where the in- 
tegration over the remaining Feynman parameters is 
implicit. 

iv) One may clearly distinguish UV and IR divergences 
in the one-loop diagrams which are contained in M s , 
see sections [| 0. All divergences in M R at d — > 4 are 
declared to be ultraviolet. At the end, all infrared di- 
vergences cancel in T^n ■ Ultraviolet divergences of the 
one-loop diagrams cancel against the divergent parts 
of the low-energy constants, as a result of which the 
quantity n is IR and UV finite. 

The calculations may be simplified considerably if one 
makes use of the fact that power counting is now pre- 
served: Nucleon (meson) propagators count as p~ x (p~ 2 )- 
The numerators can then be simplified by dropping terms 
that generate contributions beyond the order p 3 consid- 
ered here. The procedure is described in p3]. Here, one 
may use in addition that we work at threshold, where the 
pion momenta are proportional to the nucleon momenta, 



P 



qt+ = p^M^jnip. Examples of such simpli- 
fications are given in section |§|, where the triangle diagram 
is evaluated. 

We add a comment on the quark mass difference A = 
rrid — m u that occurs in the isospin violating part of 
the threshold amplitude. It shows up in two ways: first, 
through the neutron and proton masses in Feynman dia- 
grams, and second, through vertexes in the effective La- 
grangian. Consider first the latter contributions. From Eq. 
(4.1), it is explicitly seen that A does not occur at order 
p 3 , once the physical masses are used in the tree diagram 
Fig.||b. Further, in Eq. (4.1), we may replace the neutron 
mass by the proton mass at order p 3 - A still does not 
occur. Next consider loop diagrams, which are evaluated 
with the nucleon masses in the chiral limit, m p = m n = to, 
and are of order p 3 in the framework used here. Therefore, 
A cannot contribute. [Note that there is a difference be- 
tween the nucleon and pion masses in this respect: The 
electromagnetic contribution to the charged pion mass is 
of the same order in the chiral counting as the leading 
term. Therefore, it matters what pion is runni ng i n the 
loop.] We conclude that the amplitude T m in (1.2) is of 
order M% and therefore beyond the accuracy of the calcu- 
lation considered here. 



M 1 = / dxF = M s - M R 



M 



dxF , M R = 



8 The triangle diagram 

In this section we investigate the infrared-divergent 
Coulomb phase, as well as the singular terms in the real 
part of the amplitude that behave like |p| _1 in the vicinity 
of the physical threshold. At one loop, these factors arise 
when one calculates radiative corrections to the strong 
tree-level amplitudes which are obtained from the nu- 
cleon Lagrangian . For illustration, we consider the 

/°° strong diagram emerging from the Weinberg- Tomozawa 

dxF . (7.2) (WT) vertex (Fig. |a). The scattering amplitude corre- 
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sponding to this diagram is 



-, n 4+4' , n 



.1) 



The Coulomb singularities arise when one considers the 
virtual photon corrections to this vertex with a particular 
topology (vi), depicted in Fig. 0. There are two diagrams 
of this type, with the photon emitted and absorbed either 
in the initial or in the final state. The contributions of 
these diagrams to the amplitude at threshold are the same, 
so it suffices to consider one of them, and multiply the 
result by 2. The total contribution of these diagrams is 



T — 



d d k u{ P '){4+ 4'- 



2F 2 J {2ir) d i 



(to- 



m 4- jt)u(p) 



[to 2 - (p + k) 2 ] [M* - (q - k) 2 ] ' (8 ' 2) 
The numerator can be rewritten in the following manner, 

u(p')(4+ 4'- W(m+ i>+ m 4- W«(p) = 

u{p')\{2(4+ {)- ft) [ m 2 -(p + k) 2 ] 



(8.3) 



According to the discussion in section ^, we drop here 
the last term that does not contribute to the amplitude 
at 0(e 2 p). Then, after some algebra, one arrives at the 
following expression for the sum of contributions of tree 
(WT) and one-loop ( vi) diagrams to the invariant ampli- 
tude D given by Eq. dO), 



-4/(m+ jO 4 + ^4' P4 Mp) 



D 



WT 



8mF 2 



Ji 
2M1 



+ 4J 2 +4(s- to 2 )J 7 (s) 



(8.4) 



where the ellipsis stands for higher-order terms in the chi- 
ral e xpan sion. The scalar integrals that enter the expres- 
sion 



Ji 



4), are given by 
d d k 



1 



M 2 



(2ir) d i M 2 - k 2 16tt 2 
d d k 1 



uv ! 



(27r) d i [-k 2 ] [M 2 -(q- k) 2 } 
1 



J 7 (s) 



16tt 2 
d d k 



1 



(2w) d i [-k 2 ][m 2 - (p+ k) 2 ] [Ml - (q - k) 2 } ' 

(8.5) 



The representation (8.4) turns out to be very convenient 
to study the threshold behavior of the amplitude. The rea- 
son for this is that the integral J 7 (s) enters here with a 



coefficient which counts at 0(e 2 p 2 ) - so, the regular part 
of J 7 (s), which is a polynomial in the external momenta 
and masses, cannot contribute at 0(e 2 p) to the amplitude 
D. Consequently, one may use the standard dimensional 
regularization instead of infrared regularization for calcu- 
lating this integral. 

Introducing Feynman parameters, the momentum in- 
tegration in J 7 (s) can be explicitly done, 

J-y(s) — ^ _ w/ 9^ / da\da2das 



(4n) d / 2 J 
dfl-^aA [.7(a)] 



ld/2-3 



J (a) = («i + a2)(aiTO 2 + a^M 2 ) — aia 2 s . (8.6) 

One may single out the infrared divergence at d = 4 by 
rescaling the Feynman parameters according to ai,2 — * 
(1 — a^)ai 2 - the integration over the parameter a 3 then 
factorizes, and we arrive at the following expression, 



J 7 (s) 



R = a x 



1 



32tt 2 to 2 
+ (l-ai)r 2 



in 



1 da,\ 
~R~ 



1 don In R 
R 



ai(l - ax) 8 . 



•7) 



Using the formulae (C2), ( |C4|) from appendix |(], and ex- 
panding near threshold, for the sum Dy/ T + D vi we finally 
obtain 



D 



WT 



D * = : 5j£V -- > '""< 1 pi» 



+ ^ji +0 (|p|)| 
irr J 



4 = 



T - (7C UV + 8C m + 301nr + 16) + 0(r 2 ) 



rn 



where the Coul omb phase #c(|p|) is given in Eq. 

From Eq. (S.8) it is evident, that the combination 
exp -2ia0 c (| P |) (p WT + j cf with Eq (P^ j docs not 

contain the Coulomb phase at O(o). Further, dropping 
the term proportional to |p| _1 , one may read off the cor- 
responding contribution I v to the threshold amplitude 
n , see table ||, entry (i>i). Finally, we note that in the 
amplitude there is no term that behaves like In |p| near 
threshold - the coefficient Bi in Eq. ( 2.10| ) first appears at 
two-loop order |llfl. 



9 The threshold amplitude at order p 3 

We do not provide any further details of the calculations 
that can be done along the lines illustrated in previous 
sections. In this section, we combine the various pieces 
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and present the final expression for the threshold scat- 
tering amplitude T^jy. This quantity is given by a sum 
of several c ontr ibutions: to the contribution of Born dia- 
grams Eq. (t4.l|) , multiplied by the pion and nucleon wave 
function renormalization factors (Eqs. (5.7) and (6.12), re- 



spectively), one has to add the contributions from vector 
(Fig. |9|), axial (Fig. [l(]) and strong (Fig. |ll|) loops, as w ell 
as the counterterm contribution T ct displayed in (4.2), 



% N = Z„Z N T B + T V +T A + T S +T C 



(9.1) 



where 



T 



8tt 2 F 2 



6 2 2 9 

ri rr-K e 9A m \- rt 



7T 2 F 4 



(9.2) 



i=i 



The quantities I v , P A , Iq and I*, that correspond to the 
contributions from individual vector, axial and strong di- 
agrams, are listed in tables |[ and g in appendix Not 
all of the diagrams that are shown in Figs. ||, |o| and [□], 
do contribute to the threshold amplitude at 0(p 3 ). The 
diagram (s 24 ) in Fig. |Il] e.g. vanishes trivially after the 
momentum integration. Moreover, there are two reasons 
for which some of the diagrams that are formally 0(p 3 ) in 
chiral counting, start to contribute at higher order: 

i) In the one-particle reducible diagrams Fig. [l(| (ai), 
(a 2 ), (a 5 ), and Fig. [ll} (si), (s 2 ), (s 5 ), (s e ), (s 7 ), (s 8 ), 
(sg), (S17), the neutron propagator is followed by 
4-nlbu{p) (initial state), or preceded by u(p') 4irl5 (fi- 
nal state), where q^ stands either for q or for q'. For- 
mally, the combination [^75 x propagator] is of chiral 
order zero. However, putting 75 through the nucleon 
propagator changes the sign of the nucleon momen- 
tum, and the whole expression starts to contribute at 

o(p). 

ii) Doing the simplification of numerators in a manner de- 
scribed in sections ^ and || for the diagrams (03), (as), 
(ar) and (S3), (S4), (sio), (sn), it is easy to observe 
that the leading-order contributions in the numerators 
vanish at threshold. 

For this reason, in tables |(| and || we do not display 
the contributions from the above-mentioned diagrams. 

Adding all pieces together in the amplitude, we have 
checked that all UV divergences, as expected, cancel with 
the divergent parts of the LECs, whereas the IR diver- 
gences cancel among themselves. The final results is thus 
given in terms of renormalized LECs l\, fcj", d\, g\. Further, 
we split the threshold amplitude in its isosp in-c onserving 
and isospin-violating parts according to Eq. (^) . The ex- 
pressions for these parts are given by (for convenience, we 
use here the physical proton mass m p and the physical 
pion decay constant F n ) 



T° 



2FI 



9\Ml 
4m p F 2 



2_Ml 



(-2ci + c 2 + c 3 ) 



+ 



9\M* 
8m 2 p F 2 



Ml 



+ 



Ml 



16tt 2 F* 
4Ml 



3irg 



2 In- 



Mjr 



+ = K + 4 + 4 + 2dg) + 0(p 4 ) , (9.3) 



ST = 6T 2 + 5% + 0(p A ) . 



6T 3 = ST stI + ST cni + ST ct , 



ST, 



L 3 

4A 



(9.4) 



MkAk f o 33irg\ ( n , Mk 



3 327r 2 Fi 



2 In- 



;)Tr = -iMlA ( 2 + n + 8 In 2 + 12 In ^ 



32tt 2 F 2 



t _ %M^A^ 

01 3 - p2 d 5 



+ 2e'M 7! [gl+ g l-—{¥ 1 +¥ 2 )), 



(9.5) 



where ST2, ST3 are of order p 2 and p 3 , respectively. From 
the above equations one may easily check that the ampli- 
tude does not depend on the scale p. Note that here we 
have replaced the pion decay constant in the chiral limit 
F by the physical decay constant F„, 



This is the reason that l\ occurs in Eq. (3.3) 



Equations (3.4) and (9.5) represent our main result. 
This is a complete calculation of the isospin-breaking part 
of the threshold scattering amplitude at 0(p 3 ) in rela- 
tivistic baryon ChPT, by using the infrared regularization 
method of Refs. ||,||. The isospin-conserving part of 
the amplitude, given in Eq. (9.2), agrees with the one dis- 
played in appendix [A| of Ref. ]33[. A calculation of isospin 
breaking effects at order p 3 has been performed by the 
same authors in fTij] . Since they do not provide a com- 
plete analytic expression of the scattering amplitude, it is 
not possible to compare the results. 



10 The size of the low-energy constants 

In this section, we discuss the size of the low-energy con- 
stants that occur in the isospin breaking amplitude ST. 
The amplitude 5T 2 (ST 3 ) contains 4/i + f%, c\ (k\ + k r 2l 
<?g + <7g, dg). We start the discussion with ST and con- 
sider the constant c\ . It has been determined from thresh- 
old data (34| in Ref. |32|]. The uncertainty in c\ may be 
obtained from equation (20.2) in J32|, using Holder's val- 
ues jl4| for the uncertainties in the threshold parameters 



J. Gasser et al.: Ground-state energy of pionic hydrogen to one loop 



13 



and in the coupling constant g^N- We findg 

ci = -(0.93±0.07)GeV _1 . (10.1) 

For a comparison with other determinations of ci , see ta- 
ble 1 in Ref. [§§. 

Next, consider fi and f 2 that occur in the order e 2 
pion-nucleon Lagrangian. Up to terms of order 7r 4 , that 
Lagrangian is 

C^ 2) = F 2 §OW , 

= y[(/i+/ 3 )-l2 + /2T 3 ] 



In the absence of precise experimental information on 
fx, we can i) rely on order-of-magnitude estimates, or ii) 
consider model calculations. As to order-of-magnitude es- 
timates, we follow Fettes and Meifiner [B8| and write 



F 2 e'\h\ 



2tt 



l/i 



1.4GeV~ 



(10.6) 



2/ie 2 
F 2 



Q 4F 2 \V2tt°tt- 2tt+tt- 



tt+tT ■ 1 2 + ef 2 Q + 
-27T+7T- \/27r 7r + 



1 2 =diag(l,l), 



(10.2) 



where t% denotes the Pauli-matrix. From this decompo- 
sition, it is seen that fi occurs in the chiral expansion 
of the nucleon mass and in elastic pion-nucleon scatter- 
ing 7r ± p(n) — > 7r ± p(n). The electromagnetic part of the 
proton-neutron mass difference is given by the constant 
f 2 at leading order in the chiral expansion, 



2 F 2 f 2 = {m p - m n ) c 



(10.3) 



Here, we disagree with the result Eq. (12) of Ref. Q by a 
factor of 2. Numerically, we use (m p — m„) cm = (0.76±0.3) 
MeV f$7||, or 



because /i is due to a genuine photon loop at the quark 
level (we divide by 2ir rather than by An |38| to be on 
the conservative side). This estimate also confirms the 
expe ctatio n Jl3| that |/i| has the same size as I/2I, see 
Eq. ( 10.4 ). As to model calculations, we refer the reader 
to Ref. ]39|, where c\ and /i,2,3 have been determined in 
a quark model, with the result 

ex = -1.2GCV" 1 , 

F 2 (fiJ 2 J 3 ) = (-19.5 ± 1.6, -8.7 ± 0.7, 18 ±1.5) MeV . 

(10.7) 

Finally, we come to the low-energy constants that occur 
in the next -to-leading order amplitude 5%. The contribu- 
tions of the relevant LECs to S e is suppressed by one power 
of the pion mass with respect to the ones from c\ , fi and 
f 2 that occur at leading order - therefore, we expect their 
effect to be substantially smaller, because they represent 
polynomial parts of the amplitudes, not chiral singular 
pieces. The value for efij is given in table 1 of Ref. [^2| 
(we have chosen the entry with the largest uncertainty 
assigned) 



h = -(0.97 ± 0.38) GeW 



(10.4) l6M 2 d r 5 (n) = 0.04 ± 0.06 - 



2M 2 



lGeV 



We are now left with the determination of f\. The sum 



(10.8) 



m n contains the combination e 2 (f 1 + / 3 ) - the con- H ere, we have translated the dl from Ref. J?! into the 



stants /1 and fs can therefore not be disentangled from 
information on the nucleon masses. We may consider m p + 
Ttiyi as a quantity that fixes once f± is known. There- 
fore, elastic pion-nucleon scattering is the only realistically 
accessible source of information on f±. In principle, one 
may consider combinations of amplitudes that vanish in 
the isospin symmetry limit, and determine /1 from those. 
The combination 



present scheme according to Eq. (B14). Below, we will use 
Z^(lGeV) = 2.9 • 10~ 3 §. For k\ + k r 2l we invoke standard 
dimensional arguments for an estimate of its size, \k\ \ ~ 
— ^p-, and apply the same rule to <?g 



3s: 



\k\ + k 2 



16tt 2 



\9e 



16tt 2 F 2 



(10.9) 



X = T 71 



rpTT 



2T* 



(10.5) 



has this property. The tree graphs of X start at order p 2 
and contain /1 - that one may try to determine hence from 
here. Of course, one is faced with a problem of accuracy: in 
order to determine X , one needs to consider the difference 
of two large numbers, quite aside from the fact that the 
cross section n p — > ir°p is not known experimentally. It 
remains to be seen whether a combination of experimental 
data and lattice calculations could resolve the problem 
also in practice. 



[In Ref. p7| , g 7 e + g§ was estimated from an analysis of 
pion-nucleon scattering at low energies. However, in that 
reference, renormalization was treated in a manner differ- 
ent from the framework used here, and it is not clear to 
us how we can compare the couplings determined there 
with the ones used here, see appendix [b|.] As a check of 
this procedure, we apply the same esti mate to namely 
\dl I ^ 167 )i F i ■ Comparing with Eq. 
this guess is a rather generous one. 



(10. 



we see that 



6 We use for M n , m p , . . . in the following the values quoted 
in I35), in particular F„ = 92.4 MeV, |#Aj = 1.267. 



11 Ground-state energy-level shift 

We are now ready to provide numerical values for the cor- 
rection S e in the strong energy-level shift S e according to 
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Table 1. Individual contributions to 8 t . K denotes the bound- 
state correction in Eq. (2.4). See text for details. 



source 


8 t x 10 2 


|/i=o 


-5.5 


A 


±2.8 




H = 500 MeV 


At = 1 GeV 


5T£ tr 


-3.5 


-3.4 


<5T 3 cm 


0.4 


1.1 




-0.3 ±0.3 


-0.2 ±0.3 


ki ± fcj 


±0.2 


36+38 


±0.4 


K 


0.66 


vac. pol. @ 


0.48 



Eqs. (1.2) and ( |2.4| ). In addition to the LECs, we need the 
value for the scattering lengths . We take from Ref. j| 



11.1 The correction S f 



0.0883M: 



(11.1) 



The contribution from the leading term ST 2 was already 
determined in Ref. fl3|| . This contribution contains the 
three LECs fa, fa and c\. For fi, the estimate |/i| < I/2I 
was used in flgj] . The corresponding contr ibuti on to <5 e is 
±1.9 ■ 10~ 2 . Here, we use the estimate ( 10T6|) , which is 
slightly more conservative, and display in table [lj the con- 
tribution due to <572 |/ 1= o ; /1 and to the bound state cor- 
rection K. We will use below this information to estimate 
the uncertainty in the final result. 

The next order in the chiral expansion of S e is obtained 
from the second term in (2.4), evaluated with ST = ST3. 
We split this term further into contributions that stem 
from the ph oton loop, strong loops and counterterms, as 
indicated in (9A). Note that the individual terms are scale 
dependent, whereas the full result is not, of course. In 
the table, we present the values obtained from the strong 
loops and photon loops at two values of the scale, \i = 
500 MeV and /i = 1 GeV. The contributions from the 
strong loops are substantial and of negative sign. This 
large contribution is mainly due to the three graphs Fig. [ll] 
(sw), (S21) and (S22), which add 



<5T 3 str 



(Sl9),(s2l),(«22) 



1287r2.Fi 



or S e = —3.5 • 10~ 2 . We find it amusing to see that it 
is the triangle graph that generates these large contribu- 
tions. [In the self-energy type graph (sig), one has to ex- 
pand the neutral pion mass in the propagator around the 
charged pion mass. This expansion generates a triangle- 
type graph.] A graph with triangle-topology was found to 
generate a large contribution to the photo production of 
neutral pions a long time ago |^(J. In addition, that con- 
tribution is responsible for the breakdown of a so-called 
low-energy theorem in that process. Further, this graph 
also plays an important role in the low-energy analysis of 
pion-nucleon scattering |l8|,[32||. 

One might worry that the next-to-leading order cor- 
rection could be large as compared to the leading order, 
and jeopardize the chiral expansion also in the isospin 
breaking sector. However, as we already noticed, the rea- 
son for this large correction is well understood - it is due 
to the triangle graphs. We do therefore not consider this 
to be a problem. 

Finally, we display in the table the contributions from 
the counterterms in 5T ct , using the estimates displayed 
above. It is seen that the effect of the LECs at this order is 
indeed suppressed with respect to the leading order ones. 
We then add the central values at the scale /i = M p , with 
fx = g\ = k\ = 0. For the uncertainty, we add in quadra- 
ture the contributions from fx, 6 fa, Sex, k{ ± k\ , g§ ± c/g 



and from Sd r 5 . In this manner, we obtain 



-7.2 ±2.9) • 10' 



(11.3) 



This is our final result for the correction <5 e . The un- 
certainty in ( ll.Sj ) is dominated by the largely unknown 
coupling fa, as is seen from the third row in table 0. It 
does not take into account contributions from order jr in 
the chiral expansion of the isospin breaking amplitude ST. 

We note that our numerical results for the isospin- 
breaking corrections in the am plit ude cannot be directly 
confronted with those of Ref. |17[| , which refer to differ- 
ent physical quantities. Numerical values for the isospin- 
breaking effects given in Ref. ITT]] are typically of order 
~ 1 - 2 %. 



11.2 Comparison with model calculations 

Here we compare our result with the potential model cal- 
culation performed in Ref. |J , and with the evaluation |j9| 
of the couplings at leading order in a quark model, Eq. 
( 10.7 ). We start with the latter. The central values in 
(10.7) leadQ - adding the bound state correction K - to 
<5 £ = — 2-10~ 2 at order p 2 . Taken at face value, 

in our calculation: fx = 



at order p . 
this determines the coupling fa 
— 1.43GeV , in agreement with the estimate ( 10.6 ). This 
leads finally to S e — —4.3 • 10~ 2 at order p 3 . 



(11.2) 



7 In Ref. |g9j slightly different values for the pion masses 
were used. The result quoted there for the correction to the en- 
ergy shift corresponds to those values. We thank V. Lyubovit- 
skij for correspondence. 
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We now come to the potential model discussed in 
Ref. p], which predicts 



S e = (-2.1 ±0.5) ■ 10" 
We have the following comments. 



(11.4) 



The result ( 11.4 ) looks considerably more precise than 
the effective field theory evaluation. There are no un- 
known LECs occurring in the framework used in Ref. 

1- 

The leading terms in the effective field the ory f rame- 



11.3 



work - that lead to the large uncertainty in ( 
on the other hand not all incorporated [jl3|| in flll.4 ) 
An example of a QCD diagram, whose contribution 
is completely omitted, is shown in Fig. |8[ In the po- 
tential model H , it can emerge neither from the nj 
channel (there is no intermediate state with a neutron 
in Fig. H), nor from the Coulomb rescattering (the pho- 
ton in Fig. U is attached only to the proton). At the 
level of the effective Lagrangians, such contributions 
are encoded in the low-energy constants /i and fa. 
The potential model used in Ref. J(| takes into account 
terms at order p 3 , p 4 , • • • in the language of the low- 
energy expansion. We are not aware of a proof that it 
includes all of them. On the other hand, one can con- 
struct potentials that certainly do |l2| . Unfortunately, 
these potentials do not determine the LECs, because 
the LECs are used instead to pin down the potentials. 
The observation that potential models do not, in gen- 
eral, include all effects of QCD+QED is not new. Aside 
from the hadronic atom case [|[H,|l^,|llJ, we mention 
here the investigation of Fettcs and Mcifincr in ttN 
scattering 17 1. These authors pointed out that graphs 
like the one displayed in Fig. ^ (vg) generate a very 
large effect, that has not been fully accounted for in 
existing phase shift analyses. 

We conclude that the uncertainty in (11.4) is under- 
estimated - it does not reflect the systematic errors 
inherent in the method. 




Fig. 8. An example of a QCD diagram, which is not taken into 
account in the potential model of Ref. |j] . Wiggly lines denote 
gluons. 



12 Summary and conclusions 



The aim of the present work was to evaluate the 
ground-state energy of pionic hydrogen in the frame- 
work of QCD+QED. We have performed the calcula- 
tion at next -to-leading order in the low-energy expan- 
sion, relying on the method of effective field theories. 
According to the investigations of Ref. |ll| , all that is 
needed for this purpose is an evaluation of the elastic 
7r N scattering amplitude at order p 3 . 
In order to achieve this goal, we have invoked rela- 
tivistic baryon ChPT in manifestly Lorentz invariant 
form p8| , and generalized the procedure to allow for 
virtual photons. To have a powerful check on our cal- 
culation, we have in addition performed a heat-kernel 
evaluation of the ultraviolet divergences at order p 3 , 
using a consistent set of low-energy meson-meson and 
meson-nucleon effective Lagrangians that are used in 
the calculation of the 7rA?" scattering amplitudes at 

o( P 3 ). 

We have then calculated the ir p elastic scattering am- 
plitude at threshold - in the presence of isospin break- 
ing - at 0(p 3 ) in the low-energy expansion. The result 
contains unitarity corrections and counterterm contri- 
butions generated by LECs. The quark mass difference 
md — m u only enters at order p 4 , which is not consid- 
ered here. 

At leading order, the contribution to the energy shift 
is generated 
Whereas the 



L 



by three counterterms /i , /2 and c\ . 
after two can be determined from other 
sources, the size of j\ needs yet to be determined in a 
model-independent manner. 

The loop contributions, which contain mass splitting 
effects from diagrams with strong loops, turn out to be 
sizeable. Graphs with a particular topology (triangle 
graphs) turn out to be particularly important - their 
contribution to the energy shift indeed is large and 
negative, a sizeable fraction of the leading order term. 
Graphs with the same topology play an important role 
e.g. in photo production of neutral pions Efl| . 
The LECs at next-to-leading order are suppressed by 
one power of M w and thus expected to have a small 
effect on the energy shift. This expectation turns out 
to be correct for one of the couplings that has recently 
been determined [32] ] in a comprehensive analysis of 
ttN scattering. Estimating the size of the remaining 
terms with dimensional arguments, it indeed turns out 
that their contribution is about an order of magnitude 
smaller than th e lea ding order result. The final result 
is given in Eq. ( 11.3| ) 

Whereas it is true that we cannot yet provide a reliable 
error estimate of the final result, because /i is not 
known with sufficient accuracy, it is also tr ue t hat the 
calculation is systematic: The result (|9.3D-(9.5) for the 



threshold amplitude allows one to evaluate the energy- 
level shift (at this order in the low-energy expansion) 
from first principles, once one has worked out a reliable 
estimate for the LECs. 

A precise determination of the scattering lengths a+ + + 
aZ + from a precise measurement of the ground-state 
energy-level shift of pionic hydrogen has to await a 
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more precise determination of /i in our opinion. This 
fact is hidden in the potential model calculation Q, 
that quotes a very small uncertainty. As we outlined 
above, this result does not reflect all systematic uncer- 
tainties hidden in this approach: potentials in general 
do not incorporate the constraints from QCD+QED, 
unless one imposes these constraints on them Jl2| . A 
method different from effective field theories to per- 
form this matching is not in sight. We conclude that 
one is bound to know the LECs, quite independently 
of the framework used. 



We thank T. Becher, G. Colangelo, G. Ecker, T. Ericson, B. 
Kubis, H. Leutwyler, V.E. Lyubovitskij, H. Neufeld, U.-G. 
Meifiner, and M.E. Sainio for interesting discussions, and G. 
Ecker, U-G. Meifiner and H. Neufeld for useful comments on 
the manuscript. The work was fulfilled while M.A.I. , E.L. and 
M.M. visited the University of Bern. This work was supported 
in part by the Swiss National Science Foundation, and by 
TMR, BBW-Contract No. 97.0131 and EC-Contract No. 
ERBFMRX - CT980169 (EURODA<5NE), and by SCOPES 
Project No. 7UZPJ65677. M.A.I, appreciates the partial sup- 
port by the Russian Fund of Basic Research under Grant No. 
01-02-17200. 



A Meson-meson and meson-nucleon 
Lagrangians 

In this appendix, we give the operator basis and the diver- 
gent parts of the LECs of the meson-meson and meson- 
nuclcon effective chiral Lagrangians, that were used in the 
calculations. We set the charge matrices to their physical, 
constant values. 



Q= 3diag(2, 



-1) 



Q = ediag(l,0), (Al) 



and use the following notation, 

[D li ,x\ = d ll x + [r ll ,x}, 

W{0 + h.c.) W = § O W + h.c. , 
V i{0 - h.c.) V = i(<PO& - h.c.) . 



(A2) 



where X denotes any bosonic operator, and & is a fermio- 
nic field. We use the notation of Bjorken-Drell |2l| for the 
Dirac matrices. 



Table 2. Operator basis and the divergent parts of the low- 
energy couplings in the 0(p 4 ) meson Lagrangian Ji^] and in 
the 0(e 2 p 2 ) meson Lagrangian [Q. The terms that do not 
contain pion fields, and terms of order e 4 are not displayed. 



i 






1 


\ (d»Wd u U) 2 
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2Z 


5 


(x t u + W x )(Q 2 ) 


l l 7 


6 


(x^u + W x )(QUQW) 
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Table 3. Operator basis in the 0(p 2 ) pion-nucleon La- 
grangian pH ], and in the 0(e 2 ) pion-nucleon Lagrangian p^ |. 



o 



l 

4m 2 



(x+> 

(u fl u u )(D^D l/ +h.c.) 

T ftv 

X+ 
— cr^F + 

8m liv 
8m " \ r tivl 
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(Q 2 + - Q 2 -) 
(Q+)Q+ 
(Q 2 + + Q 2 -) 
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Table 4. Operator basis and the divergent parts of the low- 
energy couplings in the third-order strong pion-nucleon La- 
grangian. The operator basis is the same as the one used in 
Ref. tea] for the relativistic operators. The divergent parts are 
different from Ref. [^ij , since we use the EOM to eliminate ad- 
ditional operators that arise in the HBChPT formulation, and 
work with a different basis in the mesonic sector. For compar- 
ison with Ref. pd] , see appendix [Bl 
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Table 5. Operator basis and the divergent parts of the 
low-energy couplings in the third-order electromagnetic pion- 
nucleon Lagrangian. The operator basis is the same as the one 
used in Ref. |2q] for the relativistic operators (modulo signs 
and factors oil), whereas the divergent parts are different from 
those displayed in Ref. , see appendix [b| 
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B Divergent counterterms 



It is common and useful in ChPT to identify the diver- 
gent parts of counterterms not just within a calculation 
of a particular process, but rather in full generality, using 
the background field and heat-kernel methods. One of the 
benefits of such an approach is that it provides a nontriv- 
ial check for any particular calculation - the loop infinities 
must be cancelled by the a priori known counterterms. 

In spite of the fact that the method of Ball (4lj can 
accommodate large variety of regularizations, including 
the dimensional one, it is not straightforward to extend 
it to cover also the Becher-Leutwyler infrared regulariza- 
tion 1 18 1 3lq]. One may, however, utilize the fact that diver- 
gences encountered in the infrared regularization arc the 
same as those occurring in heavy baryon ChPT [n8lp3]. 
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In HBChPT the baryon field is split into velocity- 
dependent "heavy" and "light" components^ 



N v (x) 
H v {x) 



imv-x r>+ 



///*/•• .r i ) 



P-W{x) 



P± = -(l± 7 %) 



v 2 = l 



(Bl) 



and the heavy component H v (x) is integrated out. The 
propagator of the remaining light component is a homoge- 
neous function of the residual momentum = — mv^, 
leading to a consistent chiral power counting in dimen- 
sionally regularized HBChPT. 

The heat-kernel method in HBChPT leads to a lengthy 
and cumbersome calculation. The amount of work 
required can be reduced significantly by introduction of 
the so-called super- heat-kernel method in recent work of 
Neufeld et al [29,3(1. Following closely the procedure of 



At the next step, one expands the action up to second 
order in the fluctuations £ and r\ (Euclidean analog of the 
formula (3.1) in Ref. f3C(] ). Once the explicit form of the 
fluctuations is worked out, one can use (the Euclidean ver- 
sion of) Neufeld's master equations (3.19)-(3.21), Ref. p0[ . 
At the end, the results are transformed back to Minkowski 
space. 

After some straightforward, yet tedious, algebra one 
finds the divergent counterterm Lagrangian in a "raw" 
form, 



Neufeld, we worked out the complete divergence structure 
of the one-loop generating functional corresponding to the 
Lagrangian 



r - 4- r(p 2 ) + £0 2 ) + r 

HB 7r i i -7r ' ") 



where 



C 



(p) _Tf 



HB 



(B2) 



(B3) 



( v) 

is the HB projection of C Jr , with 



The starting point of the background field machinery 
is a decomposition of fields in the Lagrangian into back- 
ground fields and fluctuations. Here one has to distinguish 
between two different elements of the chiral coset space 
SU (2) L x SU (2) R / SU (2) v , namely u L and u R , occur- 
ring in the definitions of various fields appearing in chiral 
invariant Lagrangians, 



U 



F, 



u R u' L 



u R (d M 



iRn) u R - ut (dn - iL M ) 



ul 



u R {d ll -iR ll )u R + v} L {d ll -iL ll )u L . (B4) 



Decomposing the fields into classical background fields 
and fluctuations 



£/2F 



u R = u c \e 
N = N cl + r] 



ul 



(B5) 



and switching to the Euclidean formulation of the theory, 
which is more convenient when dealing with photons, one 
first merges the bosonic fields into 7- component objects 

c, 



(B6) 



We follow the tradition of denoting the velocity by the 
same symbol u M which was used for the external vector field. 
The meaning should be always clear from the context. 



'JV7 > 



(B7) 



£raw 



167T 2 (d-4) {12 



— (u-u) +- (u^u u ) (u^u u ) 



^{[u^u v ]T^ ) + \{u-u){ X+ )\ - (BS) 



/-•raw 



F 2 



^ 16tt 2 (d - 4) 



27 + 4Z 
20 



(u-u)(Q 2 ) 



+ 



2Z{u-u)(Q\- Q 2 _) 
\((u»Q + ) 2 + {u,Q-f 



+ 2Z{(u u Q+) 2 - (u u Q\ 2 



^(x + )(S 2 )+^(x + )(Qi-Q 2 -) 
l -(u,[C^Q + ]+u,[C^Q-]) 

75 - 80Z+168Z 2 f2/ ^ 2x 2 
50 

15 + 2Z + 12Z 2 



-F 2 (Q 2 Y 

F 2 (Q\ - Q 2 _) (Q 2 ) 



/•raw 



3 + ,Z + 24Z 2 F2{Q2+ _ Qlf 



-<[£>„, «"] [Q-,Q+])}, 



1 1 nI&KSvd,*"] 



(B9) 



16tt 2 (d - 4) F 2 {8 

1 + 5g\ 



12 ^ 



[u^, [iD^,v ■ u}} 



4~ffl 
8 

4ff a - g\ 
8 



v ■ u, [iv ■ D, v ■ u] 
S ■ u (u ■ u) 
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6g A - 6g\ + g\ 
12 

-8su + g% 



(u^S ■ u) 



—v ■ u (v ■ uS ■ u 
12 v 



S ■ u ( (v ■ u) 

4s 'a - g 3 A 



S ■ u (x+) 



2 r 



1 + 5g l 
6 

4ffi + 3 5y 5 



D F^ v 



v v +g A iS fJl v l , F^ v ,v- u 



16 



iv x e XlMUp {uMp) 



ZflS^ST] (« M [iv-D,u v ]) 



1 + 5ffl 
12 



[[i£> M ,u M ] ,vu] 



9 -fv x e x ^o (F+^u p ) -3g 2 A i (v ■ D) 3 
g A v ■ D S ■ uv -15 

16 

,,4 



((u ■ u) iv ■ D + h.c.) 



+ ^4 ^( v . u )^ iv . jD + h . c .) 

+ ^-([v D,S -u]v D + h.c.) 



([S ,i ,S u ]u li u v iv- D + h.c.) 



+ g A ([S»,S»]F +liV vD + h.c.)}N, 



(BIO) 



- - (2 + 3. 9 ^ + 24Z 5 2 ) ((Q2 _ Q 2 ^ iv . D + h c ) 

+ 1 (2 + Zg\ + \2Zg\) ((Q+f iv ■ D + h.c.)} N . 

(Bll) 



The results for £™ w and £™ w 



and the result for C r § w is consistent 
recise, our £™ w is equal 



are consistent with Eq 

(A.ll) of Rcf. g~' 
with that of Ref7~|42| (to be 

to the final result of Ref. |42fl . up to the replacement 
[Dp, u^} |x_ - 1 (x_». The result for C 1 ™ agrees with 
Stcininger's expression Eq. (3.124) up to several signs 
and factors of i [Note e.g. that the expression (3.124) in 
Rcf. [^U is not hermitean, according to the definition (2.9) 
in that work.]. 

The "raw" divergent Lagrangian, obtained so far, can 
be consistently used for renormalization of any process 
up to the corresponding order. On the other hand, it is 
perhaps not the most convenient choice, since it differs 
from the standard Lagrangians used in the field. To bring 
our result to the standard form, one has to perform the 
standard procedure, namely use the basis from Rcf. [23|] for 

( A \ — 
, and eliminate some of the terms by using equations 

of motion (EOM) for the classical fields. 

The EOM are obtained from the linear part of the 

fluctuation Lagrangian, 



= -jX- - i (x- 



AiZF 2 [Q+,Q- 



- ^T a Nr a S ■ DN + h.c.) , 



iv ■ DN = -g A S ■ uN 



(B12) 



Here, we have not included an external source for the 
nucleon field, because we are only interested in the S- 
matrix elements. The proper way of using these EOM is 
to use appropriate field transformations. For the mesonic 
Lagrangian, however, the use of the lowest order EOM, i.e. 



N (8Z - g\) (Q\ - Q 2 _) S ■ u replacing the structure [D^ u^] by the RHS of the EOM, 



^ " 167r 2 (d-4) I 2 

_ ^d( 4 + 4Z-<&) (Q + ) 2 S-u 

- 2g A (1 + Z- Zg\) (Q+S ■ u) (Q+) 
+ AZg A (1 - g\) (Q+S ■ u) Q+ 

+ 2g A (l-g 2 A )(l-2Z) (Q_5 ■ «) Q_ 

+ (1 - 3g\) (Q-v ■ u) (Q+) - 4 (Q+v ■ u) Q_ 

- 2 (l - 3g A ) (Q-v ■u)Q+-2 [Q + , iv ■ c+] 

-^(1-9^) [Q-,ivc-\ 

- 2«Q+) Q+tv- D + h.c.) 



is equivalent to the systematic performance of field redef- 
initions ]43]] . For the baryonic Lagrangian, the required 
nucleon field transformations are explicitly given in p6| . 

The resulting divergent nucleon Lagrangian corre- 
sponds exactl y to the heavy baryon projections of the 
Lagrangians (3.3) and the pertinent tables in the ap- 
pendix^. For the relativistic Lagrangians at 0(p 2 ), 0(p 3 ) : 
0(e 2 p), the simple replacements 7^75 — > 2S^, iD^ — > 
mv^Ny and u^ v — > — 2z[5 M , S v ] are all what is needed. 
Note that in appendix H we do not display terms in the 

Lagrangians and £^ p , which do not contain pion 
fields. In addition, we drop all terms at 0(e 4 ). 

We wish to note that all /3-functions give n in the ta- 
ble ||, coincide with the ones from Refs. f23|,p4|. The fj- 

9 There are misprints in that equation, not present in the 
final result displayed in Eq. (3.6) of that work. 
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functions from table || are consistent with the results of 
Ref. |2q] (the latter contains a misprint in /?n, 2/3^ — » 
—2/3^). Note that in Ref. p6| , a different basis in the 
mesonic sector was used - consequently, in order to com- 
pare, one has to transform the /3-functions, see below. Af- 
ter this is done, Pi from table | are equal to Pi (i = 1 • • • 3) 
and to Pi+i(i — 4 • • ■ 23) of table 1, Ref. [g6|. The entries 
of the table ||, necessary for a coherent result, were, to best 
of our knowledge, not presented yet. As it has been men- 
tioned above, the /3-functions for the 0(e 2 p) Lagrangian, 
which were presented in Ref. p8||, and were used in the 
numerical fits in Ref. |17|], correspond to the "raw" form 



( Bll ), rather than to the Lagrangian which is brought to 
the "standard" form in the mesonic sector by using the 
EOM. 



Finally, we comment on the relation between the two 

fp 4 ) 

most commonly used choices of operator basis in the C)f 
Lagrangian. Some of the /3-functions are changed when 
one changes the basis from |2j| to (the barred quan- 
tities belong to the latter basis) 



1 + W A 



24 



Pl8 = , 



Pl9 = , 



CT 8 = 77 



(B13) 



Furthermore, some of the renormalized LECs do differ in 
the two bases, namely 

^i8 = d r 18 (n) - -^r 4 {fi), 



4F 2 

fc>( M ) = *£(/,) + \ zfm 



(B14) 



C Integrals 

In this appendix, we present some useful formulae that 
we have used in the calculations. Note that all formulae 
given here correspond to the case of standard dimensional 
regularization. 

In the diagrams with virtual photons, the follow- 
ing definite integrals over the Feynman parameter x are 
needed [here, R = x + (1 — x) r 2 — x(l — x) s is a quadratic 
polynomial in x, and s = s/m 2 > (1 + r) ], 



/ 



dx In R = —2 + X- In h x + In - — — — - / ;vrr, , 



1 -x 
2 r lnr 



1 + r 



1 — x+ 

iira + O(a ), (CI) 



dx ■ 



1 1 



In 



R sffo ( 1 — X- 1 — x+ 

O(a ) 



1 2 

r (l + r) z cr 



(C2) 



dx — = l X— In ■ 



R sao 



1-X- 



In- 



1 — x+ 

1 + r + In r 
(1 + r) 2 



i it (x- + x + ) 
2r 



(1 + r) 3 CT 



+ 0{<j ) 



(C3) 



lni? 1 f . fl-x A 
dx— = ^- 2Li 2 i 

R 8 <7q I V 1 — X- 



4tt 2 1 i2 /1-i h 

H h- In 2 

3 2 Vl-a;_ 



-2Li 2 



x + 



In 2 



+ 2 (Ins + 2 In cr ) 



+ l ln [x- ] 



2 \ 1 — X- 



2 2 2 lnr 



(l + r) 2 ao 



"(1 + r) 



where 



x± 



<7r 4 1n[(l + rH] + 
(1 + r) 2 cto 

_ s-l + r 2 ±A 1 / 2 (s,l,r 2 ) 



2.s 



(To = — = 



A 1/2 (s, 1, r 2 ) 2|p| 



s V s 

2 



A(a;, y, z) = (x - y - z) -iyz, 
ln(l - t) 



Li 2 (z) 



dt- 



t 



(C5) 



We also recall a useful representation [}f4| o f the vertex 
functions that occur in the diagrams of Fig. |ll] (si 2 ), (513). 
These are UV and IR finite, so one may put d = 4 from 
the beginning. 



l 



d 4 fc 



(i; 

[m 2 — fc 2 ] 
1 



m| - (fc + pi) 2 ] [m§ - (fc + pi + p 2 y 



(C6) 
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J 3r 



hiPi + I32P2 
1 



h = 



dx { —7— In i?i 



1 



- — ° In i? 2 t~~ In i?3 

^2 ^3 



(C7) 



D Contributions from individual Feynman 
diagrams 



In this appendix, we list the non-vanishing contributions 
to the threshold amplitude, due to the diagrams displ aye d 
i n F igs. ^ 10 and The notation used is the on e in (9.1), 
(9.2). The quantities Ctjv.ir ar e defined in Eq. (6.4), and 



r = M^/m. 



'31 



'32 



dx 



1 



(1 - a ) 



x — 1 — ao ) In h'j - ii„ 



i? 2 

x — 1 — ao ) In i?2 + ao x 



^3 



-1 - a ] In /«'; Hi) r 



1 



(1 - go) 
A 2 



ap 
^3 



-R2 

x — 1 — ) In i?2 + 



i?3 



x — 1 — ) In i?3 + x 

A3 



(C8) 



(C9) 



Table 6. Vector-type electromagnetic diagrams I v in the in- 
frared regularization, up to and including 0(r). 



Fig. | 


J V 


M 


-| (7 Cuv + 8 C* m , + 16 + 30 lnr) 


M 


§ (7Cuv + 8C m + 16 + 30 lnr) 


(«3) 


§ (Cuv — Cm) 


M 


-£(C uv + 2Ci R -4 + 6 lnr) 


(«b) 





(««) 


-5 (-3C uv + 4-6mr) 



i?l = IE 777^ + (1 — .x) 7772 — X 0- ~ X ) Pi J 

i?2 = 2; 777^ + (1 — x) mg — X(l — X) p\ , 

R3 = X 7772 + (1 ~ X ) m 3 — X (^ ~ X ) P\ > 

P3 = (Pi +P2) 2 , 

Ai = \ 1/2 (pj,plpl) -x + ai, 

A 2 = (1 - a ) A 1/2 (p?,p^,p^) • x + a 2 , 

A 3 = -a X 1/2 (pl,pl,pl) ■ x + a 3 , 

2 2,2, ,-2 2 2\ 

ai = m 2 — m 3 + p 2 + ao (m 1 — m 2 — ) 
0-2 = &3 = "^2 — m 3 — P2 
+ a (777? - 7773 +P2 -P3) I 

pI-pl + AV^f^l^D 



a 



Pi 



(CIO) 



(Cll) 



2P? 



(C12) 



Table 7. Axial-type electromagnetic diagrams F A in the in- 
frared regularization, up to and including 0(r). 



Fig.0 




(a 4 ) 


-6rC uv + 2r - 12r(ln2 + lnr) 


(as) 


-3rC UV -It- 3r(ir + 2 lnr) 


(a 9 ) 


3rC uv - 2r + 2r (tt + 2 In 2 + 3 In r) 



The integrals ^3, /31 and J32 are not changed under scaling 
of the arguments of the logarithms by any arbitrary con- 
stant value. This property allows one to replace In Ri — > 
In Ri /m 2 in the above expressions. Note that in the thresh- 
old amplitude we have A(pf ,£>2)P§) = 0. This considerably 
simplifies the calculations. 
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Table 8. Isospin-conserving strong contributions ig an d iso- 
spin-breaking strong contributions 7^ in the infrared regular- 
ization, up to and including 0(r 3 ) and O(r), respectively. 



Fig. 




7T 


(*12) 


2 

^ r 3 (3 C uv 
+2 + 6 lnr) 


2 

-^r(3C UV 
+5 + 6 lnr) 


(*13) 


2 

-§f r 3 (3 C uv 
+2 + 6 lnr) 





(«14) 


3^r 3 (C uv + 2 lnr) 





(*16) 


i5s r3 ( — 5 Cuv 
+8 - 10 lnr) 


— 128 r (3 Cuv 
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(sie) 


2±g r 3 (-5 C uv 
+8 - 10 lnr) 





(sis) 


— — 7r r 





(sig) 


_£i 3 
48 " ' 


si 7rr 

32 " ' 




~256 r 1- & °UV 

+8 - 10 lnr) 


u 


(*2l) 


§|r 3 (9C uv + 6 
+2?r + 18 lnr) 


3s 2 4 
16 " ' 


(S22) 


13 si 3 

192 " ' 


13 9 a 

— 7T r 

128 


(S23) 


-3§4 r3 (Cuv + 2 lnr) 


3k r ( c uv 
+ 1 + 2 lnr) 



10. 



11. 



12. 



13. 



14. 



15. 



16. 



17. 



18. 



19. 



20. 



21. 
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Fig. 10. The electromagnetic axial-type diagrams. 
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Fig. 11. The strong one-loop diagrams. 



